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SUMMARY 

This paper presents a method for analyzing the dynamic behavior of two- 
dimensional impacting struts and trusses of the type presently used in the construction 
of landing gear for lunar landing vehicles. The struts and trusses may contain shock 
absorbers. The method uses lumped masses to represent the system, and two- 
dimensional motion is assumed. The equations of motion are numerically integrated to 
obtain response time histories of each of the mass points, and to allow detailed analysis 
of the dynamic behavior of the system. Application of the method to a particular strut 
and truss indicates that axial elasticity is much more important than lateral elasticity 
in the dynamic behavior. The present limitation of the method is the large amount of 
computer time required for problems in which the time period of interest is large com- 
pared with the period of the highest natural frequency of the system. 

INTRODUCTION 

Current concepts for space vehicles for soft landing on lunar and planetary sur- 
faces include landing gears consisting of legs which are trusses constructed of tubular 
struts. The truss members may contain some type of shock absorber. In order to 
assure proper functioning of such landing gears, it is essential to understand the 
dynamic behavior upon impact with a landing surface. Especially important are the 
buckling behavior of the landing gear and the effects of landing- gear elasticity on the 
landing dynamics of the vehicle. 

Previous analyses of landing dynamics (for example, refs. 1 and 2) have included 
axial but not lateral vibrations of the landing gear. Drop tests of a model of an early 
version of the lunar module have shown that the landing gear undergoes lateral vibra- 
tions during impact. The effect of lateral vibrations was not clearly understood. 
Accordingly, the study reported herein was carried out. It consists of a mathematical 
analysis of the dynamic behavior of impacting trusses containing shock absorbers with 


lateral oscillations included as a possible deformation. The analysis, which is restricted 
to two-dimensional motion, represents the truss by a finite- element lumped-mass system, 
uses Euler's method of numerical integration to obtain response time histories of each 
mass point, and thus allows detailed examination of the dynamic behavior of the truss. 

The analysis is applied to a specific strut and truss configuration to serve as an example. 

SYMBOLS 

dj,d 2 distances between vehicle center of mass and masses kj and kj + 1 

EA axial stiffness 

El flexural stiffness 

F force 

I moment of inertia 

k = 

l length of beam-column element 

m mass 

M moment 

P axial load 

t time 

X,Y axes of inertial coordinate system 

x,y distance along axes 

a slope in inertial coordinate system 

6 total shortening of beam -column element 

shortening due to bending defined by equation (14) 
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62 shortening due to compression defined by equation ( 26 ) 

?7(4) lateral deflection in local coordinate system 

rii = 71 (l) 

6 slope of free end of beam-column element in local coordinate system 

1,77 axes of local coordinate system 

p radius of gyration of cross section 


Subscripts: 


C 1 

c 2 

cm 

i 

3 , 3-1 

3 , 3+1 

k l 

kj+l 

X,Y,of 

t 


mass on main strut to which secondary strut is joined 


mass at end of secondary strut which attaches to main strut 


vehicle center of mass 


jth mass point 


jth mass point referred to first local coordinate system 


jth mass point referred to second local coordinate system 


mass point at vehicle end of main strut 


mass point at vehicle end of secondary strut 


components of forces 


total 


external 


Matrix notations: 

} column matrix 

[ ] square matrix 
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Dots over a symbol denote differentiation with respect to time. Primes denote 
differentiation in the local coordinate systems with respect to £. A bar over a symbol 
denotes a nondimensional quantity. 


ANALYSIS 

The basic method of analysis is first developed for a strut. In a later section the 
method is generalized to include a truss containing shock absorbers. 

A strut is represented by a number of point masses connected by massless beam- 
column elements which may undergo both axial and lateral deformations. (See fig. 1.) 
Each of these elements is uniform, but the length and stiffness parameters may differ 
from one element to the next. Only two-dimensional motion is allowed. 

Coordinate Systems 

As shown in figure 1, three different types of coordinate systems are used in the 
analysis. The axes of the inertial coordinate system are denoted by X and Y. In 
addition to the inertial coordinate system, there are two local coordinate systems asso- 
ciated with each mass point. The local coordinate systems are attached to the strut and 
move with it. The axes of the first local coordinate system associated with the jth mass 
are denoted by £j j_j and The origin of the first local coordinate system is 

located at mass j - 1. The j_j-axis is tangent to the strut at mass j - 1 and posi- 
tive in the direction of the jth mass. The positive direction of the 7]j^_j-axis is 90° 
counterclockwise from the positive £jj_j-axis. 

The axes of the second local coordinate system associated with the jth mass are 
denoted by j+l and ’IjJ+l* The origin of the second local coordinate system is 
located at mass j + 1. The £j ? j + i-axis is tangent to the strut at mass j + 1 and posi- 
tive in the direction of the jth mass. The positive direction of the ^jj+j-axis is 90° 
counterclockwise from the positive £jj + j-axis. 

Equations of Motion 

Newton's equations of motion for the system take the following form for 
j = 1, 2, . . ., kjt 


m j*j = F x h 

j-1 + F X iji+ l + F e,Xj 

(la) 

mjyj = F Yji 

.3-1 + Fy ),H 1 + Fe > Y i 

(lb) 

Ij&j = 

F aj,j-1 + F aj,j+1 

(lc) 
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where 


mj jth mass 

Ij moment of inertia associated with jth mass 



i ,Fy m-i’ f °M-» 


forces and the moment acting on jth mass from beam 
column j,j-l 


Fx y+i’ FY j,j+i ,Fc 


forces and the moment acting on jth mass from beam 
column j,j+l 


r e,Xj’ F e,Yj 


components of external forces acting on jth mass 


Quantities with nonpositive subscripts or subscripts greater than kj do not appear. The 
procedure for computing the quantities F Xy_i> F Yy_i’ F «j,j-1’ Fx j,j+1’ Fy j,j+r 

and is given in the following section. Since the equations which are used to 

compute these quantities are nonlinear, the values are computed numerically for each 
time step. Once the values are determined, Euler’s method is used to integrate equa- 


tions (1) one time step with the assumption that Fx; , i , 

, v«rl T nftwiniM nAncfonf <^iiT*iriiT tVlo l 

r i,i+i 


FY i,H’ F “i,i-i’ Fx i,i+i’ 


Fy, u1 , and F^j j +1 remain constant during the interval. 


Force-Displacement Equations 


To determine the quantities Fv. . , , Fy. - and Fqt • • the elastic forces 
exerted on the jth mass by the beam column connecting masses j - 1 and j must be 
determined. These forces can be determined by standard procedures. Classical beam- 
column theory gives equations relating the forces and displacements of the two ends of 
the beam column. These equations are more conveniently expressed in the first local 
coordinate system. The position of and slope at the jth mass in the first local coor- 
dinate system are given by the following formulas (fig. 1): 


Also define 



j - Xj.jJcos o-j.j + - yj_i)sin o!j_i 

(2a) 

”y-i ' -f 

cj - x H )sin a yl + (yj - yj_i)cos i 

(2b) 


0 j,M " “j ' 

(2c) 



(3) 


i 
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where ijj.j is the undeformed length of the beam column connecting masses j - 1 
and j. The quantity Sjj-1 is called the total shortening of the beam column con- 
necting masses j - 1 and j. 

For notational convenience, the subscripts are dropped with the understanding that 
the quantities are in the first local coordinate system. Note that in the first local coor- 
dinate system, the beam column connecting masses j - 1 and j may be viewed as a 
cantilevered beam column (isolated section of fig. 1). In general, three types of loads 
may be needed at the free end of this element to hold it statically in its instantaneous 
deformed configuration: a lateral force F, an applied moment M, and an end load P. 
The following procedure is a well-known method of relating the forces to the 
displacements. 


The equation for the lateral deflection of the cantilevered beam column shown in 
the isolated section of figure 1 is 


d 2 77(£) _ M t 




El 


(4) 


where 


vtt) 

lateral deflection 

El 

flexural stiffness 

Mt 

total bending moment 


(5) 

(6a) 

(6b) 


Note that 

v(0 = v t = %j- 1 

where ??j is given by equation (2b). 

The boundary conditions are 

77(0) = 0 

MD = o 


Two different cases may arise: for case 1, P = 0; for case 2, P * 0. 
Case 1; P = 0.- In this case, equation (4) becomes 



( 7 ) 
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The solution of equation (7) subject to boundary conditions (6) is 


v(£) EI 


(M + FI) 


| 2 Fj£ 
T “ 6 


At the free end 


m - <?« = 


(M + FZ) 





( 8 ) 

(9a) 


or in matrix form 


e = »'6) = gj 


(M + FZ)£ 


Fl‘ 


r~\ 


r>-h 


9 


3 

?2 


F 


w 


(9b) 


( 10 ) 


Inverting equation (10) gives the relationship between the forces and displacements 


r 

F 

< 

M 

v J 


12EI 


l± 

2 



% 

< 

0 

V J 


( 11 ) 


A convenient dimensionless form of equation (11) is 



where 



( 12 ) 


(13a) 

(13b) 

(13c) 


The longitudinal displacement of the free end of the beam due to bending of the 
beam is (ref. 3): 
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•$£'*** 


( 14 ) 


The quantity is called the shortening of the beam due to bending. 
By using rj(£) from equation (8), 




' M + F iV F 

El r " 2 El 


^ 1 

d H 


M + FlYl 
El 


2 ? 3 F (M ± FZV 4 

Ell El U + 4\EI/ 5 


(15) 


Since the end load P is zero, there is no shortening due to compression. Hence, for 
the present case, 5 = 5j. 

In dimensionless form, equation (15) becomes 


»i-s 


where 


i(M + F) 2 -iF(M + F) + iF 2 




Case 2; P f 0 .- In this case, equation (4) takes the form 

j2 


The solution of equation (18) with boundary conditions (6) is 

*0 - - §0)e^ ♦ } (**”'** * 


It follows that 


M + FI - P Vi it 

— _ — 




m = v (l) = F i- Z - == + M i i == 

1 \ p p ^ e z#7Ei + e -i^p/Ei) l p p e z\/p7Ei + e -i\IWEi 


(16) 

(17) 

(18) 


(19) 


(20) 


,11 2 \ /l FB- ^n/p/EI -Z\/p7eT\ 

0 = 77 ’(i) = F --- — — +M — ; (21) 

i p p zVp7ei -iWni l p ^i z^/ei -isjp/m 


8 




In matrix form equations (20) and (21) are 




- 




l 1 . lEI e 1 


]. fijEl _ e -l \/P/EI 


1 1 


p " PVP' e Z^7Ei + e -Z\/P7Ei P P e l\l¥7El + e -l\Jp/El 

}\J¥7 m + e -i\l¥jEi 


1 _ _ 1 _ 2 

p P e lJp/El + e -l\l¥/Ei 


1 , P e 1 


p Pi ^ \/p/EI + Q -l \Jp/EI 


r -n 




> ( 22 ) 


M 




Inversion of equation (22) gives 
r 
F 


H 


M 


^ J 


1 rj e zjp7il_ e -z^/p7Ei 

Ppi e Z\/P/EI + e -Z\/P/EI 


1 


P e Z^7EI + e -Z\^7El P 
J i M - e - 1 ^ 7 * 1 


P 6 z^7ei + 6 -z\/p7ei p 


p P|P e Z^P/EI + e -ZVP/EI 


h 


> (23) 




where 


, z\/p7ei_ p -z\/p7Ei 

Di = rs\ p e 


p2 ljEf e ^/p/EI + e -Z^P/EI P 2 e Z\/P/EI + e -Z\/P/EI P 5 
The shortening due to bending is 

\2 


(24) 


* p ~ P,,i - 1 § - o 

/M + FZ - Pr?A 2 F 2 EI 


F^Z PZ 

2 p2 ’ 4 El 


(25) 
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In this case, however, 6 j is not the only contribution to the shortening of the 
beam. The end load produces an additional shortening (or lengthening) 62 , which is 
computed by the formula for compression of a bar. 

, _ PZ 
2 EA 

The total shortening (or lengthening) is then 

6 = 6j + §2 

Equations (23) and (27) determine F, M, and P if rj^, 8, and 6 are known. 
However, P is determined only implicitly; thus, an iterative numerical solution of the 
equations is required. Because of the latter point, it is convenient to consider separately 
the two possibilities which can occur: (a) P < 0; (b) P > 0. 

Case 2(a); P < 0 .- In this case equations (23) and (27) can be put in the following 
dimensionless form: 


( 26 ) 


(27) 


r ■> 


F 

< 

)= M? 


Dl 

M 

v J 


if 

1 v 2 fi 


kZ 


1 - 


sin kZ 
cos kZ 

1 

cos kZ 


cos kZ 
1 sin kZ 


kZ cos kZ 


- 1 


r 


8 


(28) 


where 


+ ——7 + -£-r|M + F + 


(kZ ) 4 (kZ ) 4 

2F 


}(kZ ) 4 \ 2kZ / 2(kZ ) 2 

[p + F + (kZ) 2 ^]si 


F + M + (kZ) 2 m 


9 2F 2 

sin^kZ sin kZ 




M + F + (kZ) 2 m 


(k 1 )' 


k 2 EI 

EA 


(cos &-!)> + 


k-\E 

VEA 


D x = k 4 (EI) 2 Dj 


1 - 


sin 2kZ 
2kZ 


(29) 

(30) 

(31) 


and Dj is defined by equation (24). 
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where k is given by equation (30). 

It remains to be determined whether equations (12) and (16), (28) and (29), or (32) 
and (33) should be used to compute F, M, and k l for an arbitrary set of values of rj^, 
6, and 5. In other words, if rjj, 6, and 6 are known, it must be determined whether 
P = 0, P < 0, or P > 0. Once this fact is known, the proper set of equations can be 
solved for F, M, and k l. To this end, consider 62 as given in equation (26). If 
62 = 0, then P = 0; if 62 > 0, then P < 0; and if 62 < 0, then P > 0. But by 
equation (27), 

6 2 = 6 - 61 

Hence, 6 = implies P = 0; 5 > implies P < 0; and 6 < 6j implies P > 0. 

Denote the value of 61 for P = 0 by 6j(0) and denote the value of for 
P 4 0 by 6j(P). By equations (11) and (15), 6j(0) depends only on 7^ and 6. For 

the purposes of this test only, let 

6j = 6j(0) (34) 

This procedure amounts to neglecting the difference between 5j(0) and 6^(P). By 
using this assumption, 


6 = 6^0) - P = 0 

(35a) 

6 > 6^0) - P < 0 

(35b) 

6 < S^O) - P > 0 

(35c) 
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where the arrow denotes an implication. There are heuristic arguments which indicate 
that the test may be rigorously valid, but no proofs have been found. 

Note that the procedure does actually determine values of F, M, and P which 
will produce the correct values of Tfy, 6, and 6. 

Once F, M, and P are determined, the forces which the beam column j, j - 1 
is exerting on the jth mass are given by -F, -M, and -P. The components of these 
forces in the inertial coordinate system are the desired quantities (fig. 1): 


F-y. . . = -P cos at* i + F sin Q!i_i 
J" 1 J 1 

(36a) 

Fy. . - = -P sin ail - F cos aj j 

(36b) 

F„. . 1 = -M 
“3,3-1 

(36c) 

The same type of procedure is used to compute . , Fy- . and F a . . - . 

3j J J?J 

In the second local coordinate system, the position and slope of the jth mass are 

■ -( x i - x j + i) cos “j+i - (yj - yj+i) sin “j+i 

(37a) 

”i,j + i = ( x i - x i+ 1 ) sin “J+i - ( y i - y j + i) cos “)+i 

(37b) 

0 j,3+l = “j " a j+l 

(37c) 

The total shortening is 


5 j,j+l = Z j,j+! ' ^3,3+1 

(38) 

where ljj+1 is the undeformed length of the beam column connecting masses j and 
j + 1. With the understanding that all the quantities refer to the second local coordinate 
system, the equations for F, M, and P are exactly as before. Hence, equations (35) 
are used to test 6 and the proper set of equations are used to obtain F, M, and P. 
The forces which the beam column j, j + 1 are exerting on the jth mass are -F, -M, 
and -P. The components of these forces in the inertial coordinate system are (fig. 1): 

F Xj, j+ l = P cos “j+l " F sin “j+l 

(39a) 

Fy^ ^ = P sin aj + i + F cos aj + i 

(39b) 

F a . . = -M 

“3,3+1 

(39c) 
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Boundary Conditions 

The boundary conditions are handled in a straightforward manner. For example, if 
one end of the strut is pinned, then the X- and Y- coordinates of that end are held constant 
while the slope of that end is allowed to vary as determined by equation (lc). In analyzing 
a buckling test of a single strut, the upper end of the strut is forced to descend with con- 
stant velocity while the slope of that end is determined from equation (lc). In drop tests, 
the foot of the strut or truss is considered to be pinned as long as it is in contact with 
the ground. The instant that the foot tries to leave the ground, the pin is released. 

Analysis of Landing Gear Truss 

A typical two-dimensional truss with shock absorbers is shown in figure 2. This 
truss is representative of a leg of a lunar landing vehicle for an impact in which the leg 
is subjected to no out-of-plane loads. The two members of the truss are joined by a pin 
and their upper ends are pinned to the main vehicle. The two members of the truss may 
be joined at the foot, or in the section of the main strut below the shock absorber, or, as 
shown in figure 2, in the section of the main strut above the shock absorber. 

The analytical model of the truss can have up to 25 masses, including one mass for 
the center of mass of the vehicle. A shock absorber can be located between any two 
adjacent masses of the main strut and a second shock absorber can be located between 
any two consecutive masses of the secondary strut. The following definitions are needed 
in connection with the analytical model: 

Cj subscript indicating point on main strut to which secondary strut is joined 

kj subscript indicating point at vehicle end of main strut 

kj+1 subscript indicating point at vehicle end of secondary strut 

C 2 subscript indicating point at end of secondary strut which attaches to main 

strut 

A typical analytical model of the truss using 11 mass points is shown in figure 3. 

The procedures discussed previously for a strut have been incorporated into a 
procedure for dynamic analysis of the two-dimensional truss. The method is first 
extended to a truss without shock absorbers; generalization to include shock absorbers 
is done subsequently. 

Consider, first, the pinned joint connecting the two struts. The equations of motion 
for the joint are as follows: 
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( m Cl + m C2)* c l = F Xc 1 ,c 1 -l + + Fe .Xci 

(40a) 

(“cl + “c^e, - F YciJ I + F Yci _ ci _ i ♦ F e Y<;i 

(40b) 

I P 'a. = F™ , + i 

C 1 C 1 "c^Cj-l “CjjCj+l 

(40c) 

x c 2 = X Cj 

(40d) 

Vc 2 = Fc, 

(40e) 

Ic 2 “c 2 = F « C2 , C ,-1 

(40f) 


Equations (40a) and (40b) imply that the position, velocity, and acceleration of the 
joint are determined by adding the forces from the three beam columns connected to the 
joint and integrating the resulting equations. Equations (40d) and (40e) express the fact 
that point C 2 on the secondary strut coincides with point Cj on the main strut insofar 
as the x and y motion is concerned. Equations (40c) and (40f) indicate the fact that 
the moments on the main and secondary struts do not couple because of the pin. 

The motion of the center of mass of the vehicle is determined by the forces and 
torques exerted by the two pinned joints connecting the struts to the main vehicle. The 
two upper ends of the struts and the vehicle center of mass are considered to be rigidly 
connected. The motion of the center of mass of the vehicle is governed by the following 
equations (fig. 4): 


(41a) 
(41b) 

(41c) 

where excluding the impacting leg 
m cm mass of vehicle 

x cm ,y cm coordinates of center of mass of vehicle 
Icm moment of inertia of vehicle 


m cm*cm = F t,X kl + F t,X kl+1 + F e,X cm 
HWcm = F t,Y kl + F t,Y kl+1 + F e,Y cm 
Wcm = F t,X kl ( y cm - y k x ) + F t,Y kl ( x ki “ x cm) + F t,X kl+1 ( y cm " Y kl+ l) 

+ Ft ’ Y k 1+ l( Xk l +1 ' Xcm ) 
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y "F 

e ,Xcm’ e,Y cm components of external force at vehicle center of mass 
F tv ,F tv components of total force (elastic plus external) at mass j 

t,Aj 

Note that equations (41) are coupled with equations (1). 

The coordinates of masses k x and kj + 1 are determined from the vehicle cen- 
ter of mass coordinates by the following equations for the vehicle orientation shown in 
figure 4: 


x kj ~ x cm + d l cos a cm 

(42a) 

yk x = y C m + d l sin “cm 

(42b) 

x k 1+ l = x cm + d 2 cos (“cm + ^cm) 

(42c) 

Yk 1+ 1 = y C m + d 2 sin (“cm + d cm) 

(4 2d) 


where 

dj distance between vehicle center of mass and mass kj 

d 2 distance between vehicle center of mass and mass kj + 1 

0 cm angle subtended at vehicle center of mass by masses k x and kj + 1 

Procedure for Incorporating Shock Absorbers 

For many applications, it is desirable to include a load- limiting shock absorber in 
a strut. A typical shock absorber of this type contains cartridges of crushable aluminum 
honeycomb. These shock absorbers are assumed to produce a constant force while 
stroking in either direction. The inclusion of such a shock absorber in the present anal- 
ysis is accomplished as follows. The shock-absorber force is computed as a reaction 
force to the end load P. If the compressive load in the section of strut containing the 
shock absorber is less than the crush load of the cartridges, the shock-absorber force 
is set equal to the compressive load, and stroking is not considered to have occurred. 

If the compressive load is computed to be greater than or equal to the crush load, the 
shock-absorber force is set equal to the crush load, stroking occurs, and the reference 
length of the shock-absorber section is readjusted to a new shortened length. Once the 
shock absorber has shortened and then started to open, a small frictional force is applied 
until the shock absorber opens out to its original length. Then, if the shock absorber 
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continues to open, the extensional shock-absorber force is applied in a manner analogous 
to the compressive shock-absorber force. 

Note that there are two end loads associated with each section of the strut. For 
example, in figure 3, consider the section of strut connecting masses 2 and 3. There 
will be an end load P2 3 acting at mass 2 as seen in the local coordinate system 
attached to mass 3; and there will be an end load Pg 2 acting at mass 3 as seen in the 
local coordinate system attached to mass 2. Stroking is considered to occur only if 
both of these forces are greater than or equal to the crush load. 

In any attempt to analyze the response of a specific system, some modification and 
generalization of these shock absorbers will undoubtedly be necessary. The present 
shock absorber was considered in order to demonstrate the capability of this type of 
analysis to include shock absorbers of the load-limiting type. 

NUMERICAL SOLUTIONS 

Solution of Force Equations 

As stated earlier, for P * 0, equations (23) and (27) must be solved numerically 
for F, M, and P. The solution can be obtained on a digital computer by using standard 
numerical techniques. However, these techniques, which require several iterations to 
achieve the desired accuracy, result in large amounts of computer time. Therefore, the 
following procedure has been used to determine the solution. This scheme determines 
values of F, M, and P which are well within the required accuracy and can reduce 
computer time by a factor of 10 or more. From equations (26) and (27), 



Hence, 

| pz | = ea|s - 6i| 

By equation (30), 

k 2 (EI)2 » EA|6 - 6,1 

M <«) 

Now 6j, which is the shortening of the beam due to bending, changes very little 
from one time increment to the next. Therefore, as a first approximation to k l , 


16 



(44) 


« * H 

where 5j(x) is the value of 6j from the previous time increment. For case 2(a), 

6j(x) is determined from equation (29); for case 2(b), 5j(x) is determined by equa- 

tion (33). The value of kZ as determined in this manner could be used in the right-hand 
side of equation (43) to compute a second approximation, and so on until the desired accu- 
racy is achieved. However, because 6j changes only slightly from one time increment 
to the next, results have shown that the first value of k Z is sufficiently accurate and no 
further iterations are necessary. Once kZ is determined, F and M are computed 
by equations (28) for case 2(a) (P < 0) or by equations (32) for case 2(b) (P > 0). Then F, 
M, and P are determined from equations (13a), (13b), and (30), where the sign of P 
must be chosen consistent with the appropriate case. The values of F, M, and P as 
thus determined do, in fact, produce the proper values of r^, 6, and 6 to a satisfactory 

degree of approximation. 


Integration Scheme 

The numerical integration of equations (1) is carried out by using Euler's method. 
This very simple scheme has been used effectively on other problems of this same type 
(ref. 1). It has the advantages of being self- starting and of requiring only one evaluation 
of the derivatives. Since the evaluation of the derivatives requires a considerable amount 
of computing time, this latter point is believed to outweigh any advantages of more sophis- 
ticated schemes, such as Runge-Kutta, which require several evaluations of the deriva- 
tives on each step. The inclusion of higher order differences in the present scheme 
(Adams method) produces a slight increase in accuracy. The equations for Euler's 
method are 

x(t + At) = x(t) + At x(t) 
x(t + At) = x(t) + At x(t) 

Equation (lc) requires special treatment for accurate numerical integration. The 
parameter Ij which appears in equation (lc) is the rotary inertia of the jth mass. 

Rotary inertia was included not because it was important, but because equations (1) were 
the simplest way to arrive at a consistent set of equations of motion. Equation (lc) could 
be completely eliminated from the analysis. However, the elimination of this equation 
would require an iterative method to determine the Q!j values. 

When compared with motion in the x- and y-directions, the a-motion has very little 
inertial lag. The a -acceleration is much larger than the accelerations in the x- and 
y-directions. Also, the a-coordinate of a mass point oscillates at a much higher fre- 
quency than the x,y-coordinates. These high-frequency oscillations require an extremely 
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small time increment for accurate numerical integration and result in unduly long com- 
puter runs. To reduce the computer time, an impulsive damper was introduced into the 
o!-motion. The action of the impulsive damper is explained in the following section. 

Impulsive Damper 

The behavior of a both with and without the impulsive damper is shown in figure 5 
for small values of the rotary inertia Ij. Without the impulsive damper, a diverges 
very rapidly because of numerical instability of the integration scheme. If a time incre- 
ment about one fiftieth as big were used, a would oscillate nearly harmonically with a 
maximum amplitude about equal to the first negative peak of the broken line. 

The solid line shows the influence of the impulsive damper. As the solution begins, 
a is most generally not in its instantaneous equilibrium position and an acceleration 
moves a toward its instantaneous equilibrium position. At this stage of the problem, 
the two curves are identical. Since a is moving toward its instantaneous equilibrium 
position, the computation is allowed to proceed normally. In the process of moving to its 
instantaneous equilibrium position, a acquires a finite velocity increment. Therefore, 
once a reaches its instantaneous equilibrium position, this velocity will cause a to 
overshoot and begin to oscillate. In order to damp out the oscillation, the a velocity 
is set equal to zero once a reaches its instantaneous equilibrium position. This 
procedure amounts to adding a ficticious impulsive moment to the system which is just 
enough to reduce the a velocity to zero. Thus, a is prevented from overshooting 
and is actually stopped in its instantaneous equilibrium position. As the solution pro- 
gresses, the instantaneous equilibrium value of a changes. When this happens, a is 
accelerated again and the entire process is repeated. As shown in figure 5, a can have 
very high accelerations even with the impulsive damper. The high accelerations indicate 
that the instantaneous equilibrium value of a is changing rapidly and a is lagging 
behind. Once a reaches its new instantaneous equilibrium position, the high accelera- 
tions are damped out. It is essential, of course, that a remains as near as possible to 
its instantaneous equilibrium position since, for small values of Ij, a reaches its new 
position almost instantaneously. For this reason, periodic checks must be made to 
assure that a is not lagging too far behind. In the cases so far considered, a has 
been found to follow along satisfactorily for sufficiently small values of Ij. As Ij 
increases in magnitude, the damping effect becomes more pronounced. An example is 
shown in figure 6. The upper plot is a time history of the deflection of the center of a 
vertical strut with shock absorber which is oscillating in the lateral direction. The 
lower plot is a time history of the a -coordinate of the lower end of the strut. 

One additional point should be mentioned. The time increment used in the numeri- 
cal integration must be many times smaller than the period of the highest natural fre- 
quency of the x- and y-motion. This statement remains true even if more sophisticated 
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numerical integration procedures are employed. Hence, if the period of the highest 
natural frequency of the truss under consideration is short compared with the time period 
of interest then a large amount of computer time may be necessary. This condition can 
occur even when the period of interest is as short as 0.1 second. For a system with ten 
masses representing the truss and one mass for the vehicle, the computer program 
requires approximately 0.075 second per iteration on an IBM 7094 computer. This fact 
appears to be the most restrictive limitation of the method. 

APPLICATIONS 
Buckling of Beam Column 

To check the validity of the method presented in this paper, a comparison was made 
to results obtained by Hoff in reference 4. The accuracy of Hoff's results has been dis- 
cussed by Sevin in reference 5. The system which was analyzed is a beam column pinned 
at both ends as shown in figure 7. The piston at the upper end of the beam column is 
forced to descend at a constant velocity of 0.256 inch per second (0.00650 m/sec). The 
case analyzed represents a very rapid loading of a slender beam column. In the work of 
Hoff and Sevin, the beam column had a small initial curvature. Since the present analysis 
does not account for initial deviations from straightness, the beam column was considered 
to be initially straight and it was given a small lateral velocity at the center. For this 
reason, the comparison of results shown in figure 8 indicate the same general type of 
behavior but not the same numerical magnitudes. A three-mass approximation of the 
beam column was used. The comparisons are presented in terms of the nondimensional 
parameters used by Hoff. 

The dashed curves in figure 8 show the results of the present analysis for an initial 
lateral velocity of 0.5 inch per second (0.0127 m/sec). Although different values of the 
initial lateral velocity produce different curves, the results shown are typical. In the ini- 
tial phase of the problem, the inertia of the beam column causes the lateral deflection of 
the center of the column to lag behind the static values. Eventually, this inertial lag is 
overcome and the final part of the solution consists of an oscillation about the static 
curve. Note that the axial load in the beam column reaches the Euler buckling load at a 
dimensionless time of 1.0. Hence, failure of an actual beam column should occur very 
close to this time. 

It should be pointed out that equations (28) are singular if the end load P is equal 
to the Euler buckling load. However, this difficulty occurs only if P is exactly at the 
critical value, or at least extremely close to it. In actual practice, the value of P was 
never close enough to the exact critical value to cause any difficulties. The equations 
were simply integrated right through the singularity. In actual physical problems the 
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singularity should not cause any difficulties because the axial load in a beam column will 
ordinarily be kept well below the buckling load. 

The next example illustrates the onset of dynamic buckling. A beam column is 
said to buckle dynamically if the amplitude of the vibrations tend to grow without limit 
(ref. 4). Behavior of this type is shown in figure 9 for three different cases. The num- 
ber of lateral oscillations which the column undergoes before buckling is dependent upon 
the velocity of the descending piston. The midpoint lateral deflections of the center of 
the beam column for three different values of the velocity of the descending piston are 
shown in figure 9. Prior to the onset of dynamic buckling, the number of oscillations 
which are predicted from the formula given by Hoff are approximately as follows: 

1/3 oscillation for a velocity of the descending piston of 0.256 inch per second 
(0.00650 m/sec); 2/3 oscillation for a velocity of the descending piston of 0.128 inch 
per second (0.00325 m/sec); 1^- oscillations for a velocity of the descending piston of 
0.064 inch per second (0.001625 m/sec). In figure 9, the horizontal dashed lines indicate 
the region of normal vibrations. As the data in figure 9 indicate, the present analysis is 
in agreement with the values predicted by Hoff. For a velocity of the descending piston 
of 0.256 in./sec (0.00650 m/sec) buckling occurs before the midpoint of the beam column 
begins to swing back toward its initial position; for a velocity of the descending piston of 
0.128 in./sec (0.00325 m/sec), the midpoint of the beam column begins to swing back and 
then buckles out; for a velocity of the descending piston of 0.064 in./sec (0.001625 m/sec), 
the midpoint swings out, comes back, and finally buckles out on the opposite side. 

Buckling of Impacting Strut 

Figure 10 shows a strut with a shock absorber and its analytical model. The strut 
is of the general type presently being considered for the landing gear of the lunar module. 
The strut is made of aluminum, is 12.3 ft (3.749 m) long and 6 inches (0.1524 m) in 
diameter, and has a wall thickness of about 0.03 inch (0.000762 m). At the upper end of 
the strut there is a heavy mass which weighs approximately 16 percent of the Euler 
buckling load. The shock absorber is a capsule of aluminum honeycomb with a crush 
load of 3800 pounds (16903.24 N). Figure 11 shows a typical load-deflection curve which 
results from the analytical model of the shock absorber as described in a previous sec- 
tion. Note that the assumption of small deflections is valid because the length of the 
shock-absorber section is readjusted to a new shortened length as the shock absorber 
crushes. 

Figure 12 shows the lateral displacements of the center of the strut without shock 
absorbers for various values of the impact velocity. The strut is dropped straight down 
and the bottom of the strut is stopped instantaneously when it touches the ground. The 
center masses of the strut are given a small lateral velocity in order to induce bending. 

At t = 0, the bottom of the strut first contacts the ground. 
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For an impact velocity of 1 ft/sec (0.3048 m/sec) the maximum axial load in the 
strut is well below the buckling load and there is no pronounced effect on the lateral 
vibrations. With an impact velocity of 2 ft/sec (0.6096 m/sec) the maximum axial load 
slightly exceeds the Euler buckling load. Before buckling begins, however, the downward 
velocity of the heavy mass at the upper end of the strut is stopped at approximately 
0.02 second. As the strut begins to rebound, the upward displacement of the heavy mass 
at the upper end of the strut begins to relieve the compressive load. For impact veloc- 
ities greater than 2 ft/sec (0.6096 m/sec), the effect on the lateral displacements is 
clearly evident in figure 12. 

Figure 13 shows the lateral vibrations of the strut with shock absorbers impacting 
at 2 ft/sec (0.6096 m/sec). The center of the strut reaches its maximum lateral dis- 
placement on the first vibration and then oscillates with decreasing amplitude. The down- 
ward motion of the large mass is not stopped until about 0.076 second. Buckling is pre- 
vented by the shock absorber, which limits the axial load on the strut. The shock 
absorber begins stroking at approximately 0.004 second and stops stroking at approxi- 
mately 0.076 second. The lateral deflection of the strut without shock absorbers 
impacting at 1 ft/sec (0.3048 m/sec) is shown by the dashed line in figure 13 for com- 
parison. For the strut without shock absorbers impacting at 1 ft/sec (0.3048 m/sec), 
the downward motion of the heavy mass is stopped at approximately 0.021 second. 

Effect of Lateral Elasticity on Impacting Strut 

As shown by the tests and supporting analysis of a 1/6-scale model of a lunar- 
module-type spacecraft described in reference 2, stored elastic energy can significantly 
affect the landing dynamics of lunar-landing vehicles. However, the analysis reported 
in reference 2 did not include lateral elasticity. 

To determine the effect of lateral elasticity, the strut described in the previous 
section was impacted at 2 ft/sec (0.6096 m/sec) with varying amounts of axial and lateral 
elasticity. The initial lateral velocities of the center of the strut were varied from case 
to case in order to produce different amounts of maximum lateral displacement of the 
center of the strut. A maximum lateral deflection of 1 inch (2.54 cm) corresponds to an 
axial stress which exceeds the manufacturer's recommendation for the ultimate stress. 

For the strut considered, the results indicate that axial elasticity is much more 
important than lateral elasticity in the dynamic behavior of the impacting strut. In fig- 
ure 14, the maximum upward momentum of the rebounding heavy mass at the upper end 
of the strut is plotted against the maximum lateral amplitude of the center of the strut 
for different values of axial and lateral stiffness. As figure 14 shows, reducing El by 
a factor of 2 or 3 has very little effect on the maximum upward momentum. However, 
reducing EA by a factor of 2 increases the maximum upward momentum by almost 
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40 percent. The figure also indicates that lateral vibrations have only a very slight 
effect on the maximum upward momentum. 

In figure 15, vertical-acceleration time histories are shown for four of the extreme 
conditions presented in figure 14. Here also it can be seen that the reduction in El 
produces only a small change in the vertical-acceleration time history whereas the 
reduction in EA produces a much larger change in the curve. 

Impact of Vehicle With Truss- Type Landing Gear 

Figure 16 shows a two-dimensional vehicle with a truss-type leg. The truss con- 
sists of two nearly equal struts joined at the foot. The struts are essentially the same 
as the strut considered above. The vehicle mass is approximately that of the lunar 
module at the moment of impact on the lunar surface. Lunar gravity is used. The 
vehicle impacts on one leg only. Calculation is stopped prior to impact of another leg. 

An impact velocity of 1 ft/sec (0.3048 m/sec) is used. This velocity is just enough to 
build up the loads in the struts to the crush loads of the shock absorbers (3800 lb or 
16 903.24 N). 

In figure 17, the vertical-acceleration time history of the center of mass of the 
vehicle is shown for three different cases. For the first case, the full values of El 
and EA are used for both struts. The second case shows the effect of reducing the 
El of both struts by a factor of two while keeping EA at its full value. As the figure 
shows, these two curves are practically identical. The third case shows the vertical- 
acceleration time history of the vehicle center of mass when the EA of both struts is 
reduced by a factor of two and El is maintained at its full value. In this case, the 
curve is considerably changed and exhibits the same characteristics as for a single strut 
dropped straight down. The horizontal-acceleration response time histories show 
effects very similar to the vertical acceleration. 

These results show that for the vehicle considered, a change in the bending stiff- 
ness by a factor of two produces essentially no change in the acceleration time history 
of the center of mass, and thus indicates that the motion of the vehicle is not substan- 
tially affected by the lateral vibrations of the landing gear. However, a change of the 
same amount in the longitudinal stiffness results in a much different acceleration time 
history of the vehicle center of mass. This result indicates that for the case analyzed, 
the longitudinal elasticity of the landing gear is much more important in the impact 
dynamics of the vehicle than the lateral elasticity. 

CONCLUDING REMARKS 

A method has been developed for analyzing the dynamic behavior of two- 
dimensional impacting struts and trusses which may contain load-limiting shock 
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absorbers. Lumped masses and finite elastic elements are used to represent the sys- 
tem and response time histories are obtained for each mass point. Application of the 
method to an impacting strut of the general type presently being considered for the 
landing gear of the lunar module indicates that the impact dynamics is essentially unaf- 
fected by variations in lateral elasticity. However, the impact dynamics is affected by 
variations in axial elasticity. Thus, axial elasticity is much more important than lateral 
elasticity in the dynamic behavior of the impacting strut. A similar conclusion is 
obtained for a vehicle with a specific truss-type landing gear constructed of these struts. 

The present limitation of the method is the large amount of computer time 
required for systems which have a high natural frequency with a period which is very 
short compared with the time period of interest. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., November 28, 1967, 
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Figure 2.- Typical two-dimensional truss with shock absorbers. 
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Figure 8.- Comparison of three-mass buckling test with results obtained by Hoff. 
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Figure 10.- Strut with shock absorber and its analytical model. 
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Figure 13.- Lateral deflection of center of impacting strut with and without shock absorbers. 
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Figure 14.- Effect of longitudinal and lateral stiffness on maximum upward momentum of vertically impacting strut with shock absorber. 
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Figure 15.- Effect of stiffness on vertical-acceleration time history of upper mass of strut impacting at 2 ft/sec (0.6096 m/sec). 









